Self-assembly of endohedral metallofullerenes: A decisive role of cooling gas and metal-carbon bonding

The endohedral metallofullerene (EMF) self-assembly process in Sc/carbon vapor in the presence and absence of an inert cooling gas (helium) is systematically investigated using quantum chemical molecular dynamics simulations. It is revealed that the presence of He atoms accelerates the formation of pentagons and hexagons and reduces the size of the self-assembled carbon cages in comparison with analogous He-free simulations. As a result, the Sc/C/He system simulations produce a larger number of successful trajectories (i.e. leading to Sc-EMFs) with more realistic cage-size distribution than simulations of the Sc/C system. The main Sc encapsulation mechanism involves nucleation of several hexagons and pentagons with Sc atoms already at the early stages of carbon vapor condensation. In such proto-cages, both Sc–C σ-bonds and coordination bonds between Sc atoms and the π-system of the carbon network are present. Sc atoms are thus rather labile and can move along the carbon network, but the overall bonding is sufficiently strong to prevent dissociation even at temperatures around 2000 kelvin. Further growth of the fullerene cage results in the encapsulation of one or two Sc atoms within the fullerene. In agreement with experimental studies, an extension of the simulations to Fe and Ti as the metal component showed that Fe-EMFs are not formed at all, whereas Ti is prone to form Ti-EMFs with small cage sizes, including Ti@C28-Td and Ti@C30-C2v(3)

Stacked Antiaromatic Porphyrins

Aromaticity is a key concept in organic chemistry. Even though this concept has already been theoretically extrapolated to three dimensions, it usually still remains restricted to planar molecules in organic chemistry textbooks. Stacking of antiaromatic π-systems has been proposed to induce three-dimensional aromaticity as a result of strong frontier orbital interactions. However, experimental evidence to support this prediction still remains elusive so far. Here we report that close stacking of antiaromatic porphyrins diminishes their inherent antiaromaticity in the solid state as well as in solution. The antiaromatic stacking furthermore allows a delocalization of the π-electrons, which enhances the two-photon absorption cross-section values of the antiaromatic porphyrins. This feature enables the dynamic switching of the non-linear optical properties by controlling the arrangement of antiaromatic π-systems on the basis of intermolecular orbital interactions

Accelerating the density-functional tight-binding method using graphical processing units

Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility (OLCF) with 1 to 6 NVIDIA Volta V100 GPUs per computer node, and (2) an in-house Intel Xeon computer with 1 to 2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2- and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters

Mapping behavioral specifications to model parameters in synthetic biology

With recent improvements of protocols for the assembly of transcriptional parts, synthetic biological devices can now more reliably be assembled according to a given design. The standardization of parts open up the way for in silico design tools that improve the construct and optimize devices with respect to given formal design specifications. The simplest such optimization is the selection of kinetic parameters and protein abundances such that the specified design constraints are robustly satisfied. In this work we address the problem of determining parameter values that fulfill specifications expressed in terms of a functional on the trajectories of a dynamical model. We solve this inverse problem by linearizing the forward operator that maps parameter sets to specifications, and then inverting it locally. This approach has two advantages over brute-force random sampling. First, the linearization approach allows us to map back intervals instead of points and second, every obtained value in the parameter region is satisfying the specifications by construction. The method is general and can hence be incorporated in a pipeline for the rational forward design of arbitrary devices in synthetic biology

Recent Developments in the General Atomic and Molecular Electronic Structure System

A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods such as the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock, and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as it is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite such as GAMESS is maintained and developed is considered. Future developments are briefly summarized